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ABSTRACT 

Temporally-resolved electron density measurements of solar flare plasmas are presented using data 
from the EUV Variability Experiment (EVE) onboard the Solar Dynamics Observatory (SDO). The 
EVE spectral range contains emission lines formed between 10^-10'' K, including transitions from 
highly ionized iron (>10 MK). Using three density-sensitive Fe XXI ratios, peak electron densities of 
j^Oii-2-ioi2 i cm~'^ were found during four X-class flares. While previous measurements of densities 
at such high temperatures were made at only one point during a flaring event, EVE now allows 
the temporal evolution of these high-temperature densities to be determined at 10 s cadence. A 
comparison with GOES data revealed that the peak of the density time profiles for each line ratio 
correlated well with that of the emission measure time profile for each of the events studied. 
Subject headings: Sun: activity — Sun: corona — Sun: flares — Sun: UV radiation 
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1. INTRODUCTION 

Solar flares are generally considered as increases in the 
X-ray luminosity on the Sun due to changes in the tem- 
perature and density of the coronal plasma. The increase 
in temperature is readily inferred from the presence of 
high-temperature (>10 MK) emission lines in solar flare 
spectra. However, precise values of the coronal elec- 
tron density (iVg) have been more difficult to ascertain. 
Quite often these densities are estimated from broadband 
continuum measurements which can yield values of the 
flare emission measure (EM = / JyN^dV). Deconvolv- 
ing density values from the emission measure requires a 
knowledge of the volume of the emitting plasma {V, esti- 
mated from imaging data) and a possible filling factor (/, 
usually assumed to be unity). Accurate measurements 
of the electron density are important for understanding 
both the heating and cooling of fiare plasmas, and deter- 
mining the mechanisms responsible. 

A more reliable derivation of electron densities can 
be made using density-sensitive line ratios, which do 
not require prior knowledge of the emitting volume nor 
the filling factor, under the assumption that one of 
the lines arises from a metastable level. While there 
have been many studies that have presented the den- 
sity struct ure of active regions at coronal temperatures 
/~l- 2 MK:IGallagher et al] [200llWaF ren fc WinebaTiia 
120031 iMilligan et al.l 120051 ) there have b een fewer of the 
coronal pla sma density du r ing fl ares. iMcKenzie et al.l 
(|1980H and IDoschek et all (|1981[) used data from the 
SOLEX instrument onboard P78-1 to derive time- 
dependent flare densities from an O VII ratio (2 MK) and 
found that values reached a maximum of ^10^^ cm~^ at 
the pe ak o f the impulsive phase. More recently, IMilliganI 
(j2011|) and lGraham et all (|2011l ) measured electron den- 
sities of 10^^-10^^ cm~'^ at flare footpoints during their 
impulsive phases from Fe XII (1.4 MK), Fe XIII (1.6 MK) 
and Fe XIV (1.8 MK) line ratios using dat a obtained from 
the EUV Imaging Spectrometer (EIS; iCulhane et al.l 
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I2007D onboard Hinode. IWarren et all ()2008f ) identified 
several Ca XV line pairs (formed at 4 MK) in the EIS 
wavelength range which ar e sensitive to densitie s in the 
range 10^-10" cm^^^ while iFeldman et al. ' (2008) identi- 
fied Ti, Cr, and Mn lines formed above 10 MK which are 
sensitive to densities from lO^^-lO^^ cm~'^. However, due 
to the telemetry and planning restrictions implemented 
on Hinode observations, fiare data from these lines are 
unlikely to become available. 

Density diagnostics of flaring plasma using emission 
lines formed above ^^10 MK have been even more elu- 
sive in recent decades. Among the first to identify 
such lines were lKastner et ahl ()1974D using data from the 
Goddard Space Flight Center (GSFC) scanning spectro- 
graph onboard the Orbiting Solar Observatory (OSO)- 
5 satellite. These transitions were primarily from Fe 
XIX-Fe XXIII and found in the 9-14.5 nm po r tion o f 
the EUV sp e ctrum . Following IDoschek etall ()1973[) . 
IMason et al.l ()1979( ) identified several Fe XXI line pairs 
in this wavelength range that could serve as density di- 
agnostics during flares for values of between 10^"'^- 
10^^ cm~^. A preliminary application to OSO-5 data 
yielded typical flare densities <10^^ cm~^ for tempera- 
tures --10 MK. IMason et aTl ()1984[ ) later refined these 
values, finding 4x10^^ cm~'^ for temperatures between 
10^'^ and 10 '^'^ K. A similar de tailed study of electron 
densities by iLawson fc Peacock (1984) using the OSO- 
5 solar fiare spectra found inconsistent values of 10^^'^, 
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and lO^'^-^ cm"'^ from Fe XX, Fe XXI, and Fe 



XXII line ratios, respectively. However, the OSO-5 in- 
strument took 2.5 minutes to scan the 9-14.5 nm wave- 
length range, meaning that each of the lines in a given 
ratio were recorded at different times, during which the 
flare was likely to have evolved. The instrument sensi- 
tivity as a function of wavelength was also uncertain. 

Density diagnostics with improved spectral and tem- 
poral resolution were later obtained from the Bent Crys- 
tal Spectrometer (BCS) onboard the Solar Maximum 
Mission. This was a soft X-ray spectrometer that ob- 
tained high-resolution spectra in small wavelength in- 
tervals near 0.19 and 0.31 nm. An anal ysis of Fe XX 
lines in BCS spectra was performed by iPhillips et al.l 
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Fig. 1. — Top panel; The 8-16 nm portion of the EVE spectrum showing the presence of high-temperature Fe lines (XVIII- XXIII) at the 
peak of the X6.9 flare that occurred on 2011 August 9 (solid line), and from a quiescent time before the flare (dotted line). Bottom panel: 
Expanded view of the 12-15 nm range which contains the four Fe XXI lines used in this study. Overplotted are the fits to each of the lines 
as well as the background level, indicated by the horizontal line. 



(|1983f l to obtain electron density values. In most 
cases values of the electron density were found to be 
~10^^ cm~'^ and were therefore indistinguishable from 
the low density limit. However, there has been one re- 
port of solar flare de nsities as high as 10"'^^ cm""^ from Fe 
XXII hue ratios by IPhillips et all ^996). In the stel- 
lar case, observations made with the Extreme Ultra- 
violet Explorer also revealed coronal electron densities 
as high as 10^'^ cm~^ during the most energetic events 
(|Monsignori Fossi et al.|[l996D . Coronal electron densi- 
ties of this magnitude, if confirmed, could have signifi- 
cant implications for fiare physics. 

Earlier studies of flare densities using high-temperature 
line ratio techniques have each focused on a single time 
interval during the fiare, often with an integration time 
of several minutes. Any time-resolved investigations were 
undertaken using lines formed at quiescent coronal tem- 
peratures. Here we present time profiles of electron den- 
sity, determined using line ratios with formation temper- 
atures in excess of 10 MK obtained using data from the 
EUV Variability Experiment (EVE: iWoods et al.|[20Tol) 
instrument onboard the Solar Dynamics Observatory 
(SDO). Section [5] gives an overview of the EVE instru- 
ment, the emission lines within the spectral range under 
consideration and a description of the methods used to 
analyse the data. In Section[3]we present results obtained 
from four X-class flares, while Section 2] summarizes the 
results and conclusions, and discusses the implications. 

2. EVE OBSERVATIONS AND DATA ANALYSIS 

EVE acquires full-disk (Sun-as-a-star) EUV spectra 
every 10 seconds over the 6.5-37 nm wavelength range us- 



ing its MEGS-A (Multiple EUV Grating Spectrograph) 
component with a near 100% duty cycle. The 9-16 nm 
portion of this wavelength range contains many emission 
lines from transitions in high-temperature (>10 MK) 
Fe ions (XVIII-XXIII), as shown in the top panel of 
Figure [U These high-temperature ir on lines dominate 
the E VE spectrum during a flare (see iChamberlin et aH 
[20ll) . 

Although the 9-16 nm range includes several species 
of high-temperature Fe lines, its coarse resolution 
(^0.1 nm) means that many of the observed line pro- 
files were blended with other emission lines, or were too 
weak to be detected above the level of continuum emis- 
sion. After investigating all possi ble density-s ensitive line 
pairs for Fe XIX-XXII listed in .Phillips et all (|2008i: page 
171, Figure 6.7), only three pairs of Fe XXI lines were 
deemed to be reliable (see bottom panel of Figure [H . 
These agree with those identified by iM ason et al.l (|1979l : 
Il984( ) and will be discussed in more detail in Section [2Tl 

The identification of the emission lines present in the 
EVE spectra was performed by visually comparing the 
central wavelength of the ob served emission fea tures with 
the CHIANTI (version 7; ILandi et al.l 120121 ) line list. 
Synthetic line profiles for this list were generated using 
the ionization equ ilibrium files o f .Brya ns et al. (2009), 
the fiar e DEM of [D ere fc Cookl (|i"979() . coronal abun- 
dances (jFeldman et al..il992i) . initial densities of 10^^ and 
10^^ cm~'^, and at the EVE spectral bin size of 0.02 nm. 



2.1. Fe XXI lines 
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The strongest Fe XXI emission line in th e wavelength 
range studied is that at 12.875 nm (see also lMason et alj 
119841) . This is unblended and common to several 
density-sensitive line pairs, specifically 12.121/12.875, 
14.214/12.875, and 14.573/12.875 (see Figure [21). The 
electron density of the flare plasma can be determined 
directly from the ratio of the peak intensities. For each 
density-sensitive line pair, the flux ratio as a function 
of Nf, at the temperature of maxi mum ionisation frac - 
tion, which is -12 MK for Fe XXI (iBrvans et al.ll2009f ). 
was calculated using CHIANTI. Lines were fitted with a 
Gaussian function and the peak, rather than integrated, 
flux used to calculate the ratio values from which the 
densities were derived. This was to avoid including any 
potential weaker line emission in the wings of the feature 
of interest to the total line flux. 

The Fe XXI hue at 14.214 nm is blended with an- 
other Fe XXI transition at 14.228 nm, and was ob- 
served as a single emission feature due to their central 
wavelengths being separated by only 0.014 nm, which is 
within one wavelength bin of the EVE spectra. There- 
fore, the density-sensitive ratio was determined from the 
sum of the flux of these two lines. This results in the 
(14. 214-1-14. 228)/12. 875 ratio being less sensitive to lower 
densities than if the lines were resolved (dotted line in 
Figured]). 

The 12.121 nm line was weak in many cases and is 
blended in the blue wing of Fe XX 12.184 nm but is in- 
volved in a ratio which is sensitive to a wide range of den- 
sities (solid line in Figure [2]). However, measurements of 
the irradiance were uncertain due to the presence of the 
stronger Fe XX line except in the largest flaring events. 

The Fe XXI line at 14.573 is blended with a lower tem- 
perature line which is present in the pre-flare spectra. By 
subtracting out a pre-flare profile from the flare spectra 
effectively eliminates any contribution from the unidenti- 
fied feature during the SXR peak of the flare, as any en- 
hancement from lower temperature lines are not thought 
to be significant until later in the flare (for coronal lines) 
or during the impulsive ph ase (for chromospheric lines; 
see ChambcrlineLaLlllOll). 

Each line ratio has a low and high density limit, beyond 
which the ratio is not suitable for determining electron 
densities. The low density limit for each of the ratios was 
taken to be 10^^ cm"'^. 

2.2. Background Subtraction 

As some of the Fe lines under investigation are rel- 
atively weak, an accurate background subtraction pro- 
cess is cruci al for determin i ng th e line flux solely due 
to the flare. iMilligan et al.l ()2012[ ) recently pointed out 
that the entire MEGS-A wavelength range includes the 
underlying free-free continuum emission (which is coro- 
nal in origin) that becomes enhanced during solar flares, 
steepening at shorter wavelengths. On top of this lies a 
pseudo-continuum formed by the close proximity of the 
many high-temperature Fe lines. However, as the tran- 
sitions here lie within 2.5 nm of each other, it is rea- 
sonable to assume a constant background across all four 
lines. Therefore, the entire 12-15 nm range was fitted 
with 17 Gaussian profiles plus a constant background at 
each time step throughout the flare (see bottom panel of 
Figure [1]). 
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Fig. 2. — Theoretical line ratios for three Fe XXI line pairs cal- 
culated using the CHIANTI atomic database. 

2.3. Uncertainties 

In the current release of the EVE data (level 2, version 
2), the uncertainties on the irradiance values have not 
yet been included and the validation efforts are ongoing. 
As such, each line proflle was flt with a Gaussian proflle 
with constant weighting. To calculate uncertainties in 
the measured densities, the errors in the line intensities 
(as derived from the least squares flt of the Gaussian pro- 
file) were added in quadrature to give upper and lower 
limits on the fiux ratios. There were then converted into 
upper and lower limits of the density for each time inter- 
val. If the lower limit was found to be less than that of 
the low density limit for a given line pair, then the lower 
limit was set equal to lO'^^ cm""^. The electron densities 
corresponding to these limits then provided the range of 
likely values for each line pair. 

3. RESULTS 

The methods described in Section [2] were applied to 
four X-class flares observed by EVE during 2011 and 
2012, including the X6.9 event that occurred on 2011 
August 9 which is the largest flare observed by EVE to 
date. In the GOES X-ray light curve, emission is ob- 
served to increase starting at approximately 07:50 UT, 
and the peak X-ray flux occurs at 08:05 UT. Electron 
densities for this flare were derived from the three Fe XXI 
line ratios described in Section 12.11 from approximately 
08:02 UT until 08:14 UT, during which the values were 
above the low density limit of the line ratios. The elec- 
tron density profiles derived from the three line pairs are 
shown in panels b-d of Figure [3] (solid lines), where the 
dotted lines represent the upper and lower limits of the 
measured density. 

Results from each of the line ratios show a similar trend 
throughout the main phase of the flare. The electron 
density increases above the low density limit of the line 
ratios at approximately 08:02 UT and reaches a peak 
at around 08:05 UT. There is then a slower decline in 
density, lasting a period of approximately 9 minutes, at 
which point the density had decreased to a value be- 
low the limit of the line ratios. Each of the three line 
pairs gave consistent peak electron density values of ap- 
proximately 10^^ cm~^. The precise values and times 
of peak Nf. for each line pair is listed in Table [T] along 
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TABLE 1 

Values and times of peak electron densities and peak emission measures for each of the flares presented in this study. 
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with their uncertainties. For comparison, the GOES data 
for each event are also plotted to show the relative tim- 
ing of the flaring emission (Figure |3^) . From the GOES 
data an emission meas ure profile can be de rived using the 
methods described in iWhite et al.l ()2005t ) . The time of 
peak electron density is in good agreement with the time 
of peak emission measure as determined independently 
from GOES observations (08:05:20 UT, also listed in Ta- 
ble [T]), as denoted by the vertical dashed line. 

This analysis was repeated for three other X-class 
flares: the X5.4 on 2012 March 7, the X2.2 on 2011 
February 15, and X2.1 on 2011 September 6. For 
all events, consistent peak density values of 10^^'^- 



Flux GOES If) 3 sec 



10 



11.9 



cm were obtained from each line pair, although 



reliable values could not be measured for the X5.4 flare 
using the 12.121/12.875 ratio (see Table [1]), which may 
have been due to the dominance of the neighbouring 
Fe XX line. Despite this limitation, time profiles of 5- 
30 minutes were obtained for each flare. 

4. CONCLUSIONS 

In this Letter we present techniques for determin- 
ing time-dependent measurements of the (coronal) elec- 
tron density using high-temperature (~12 MK) density- 
sensitive line ratios during four X-class solar flares. Pre- 
vious values of electron densities at flare temperatures 
have only been presented for a small sample of events, at 
a single time dur ing those flares with integrati on times of 
several minutes (iPhillips et al.lll98l .Mason et al.lll984 
iLawson fc Peacocklll984[ ). EVE now allows the time pro- 
flies of electron densities at flare temperatures to be de- 
termined, ultimately, for a statistically significant sample 
of events throughout the course of solar cycle 24, due to 
the -100% duty cycle of EVE. 

After investigating many density sensitive line pairs 
from lines in Fe XIX-Fe XXII, and accounting for 
blends within the EVE spectral resolution, the most re- 
liable ratios were found to be Fe XXI 12.121/12.875, 
(14.214-|-14.228)/12.875, and 14.573/12.8 7 5, in agree- 
ment with those identified bv lMason et al.l (|1984[ ). Each 
of these line pairs is sensitive to densities in the range 
10^^-10^^ cm~^. Applying these diagnostics techniques 
to four X-class flares revealed peak densit ies of 10^^'^- 
10^^'^ cm~'^, in broad agreement with iMason et al.l 
11984). 

The highest densities measured were for an X6.9 flare 
and were of the order of 10^^ cm^"^. The fact that consis- 
tent values were obtained using each of the line pairs indi- 
cate that any potential line blending did not signiflcantly 
affect the measurements. While comparable values were 
obtained for the other X-class flares analyzed, no values 
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Fig. 3. — Electron density profiles from three Fe XXI line ratios 
during the X6.9 class flare on 2011 August 9. Panel a shows the 
GOES lightcurves for the event in the 0.05-0.4 and 0.1-0.8 nm 
channels. Panels b-d show the time profiles of electron density 
determined using the 14.573/12.875, (14.214-|-14.228)/12.875, and 
12.121/12.875 Fe XXI ratios, respectively (solid lines). Dotted lines 
denote the upper and lower density limits. The vertical dashed 
line in each panel indicates the time of peak emission measure as 
derived from GOES data. 



of electron density above the low density limit could be 
determined for any M-class flares that were investigated. 
This limits the suitability of the line ratios investigated to 
X-class events where the electron density is greater than 
10^"'^ cm""^, except perhaps in extreme cases. These re- 
sults tentatively suggest that electron densities are larger 
in solar flares with a higher peak X-ray flux. 

Perhaps more significantly, from the time profiles of 
the density evolution, the time of maximum density was 
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found to coincide with that of the peak emission measure 
as determined from GOES data. The GOES 0.05-0.4 
and 0.1-0.8 nm passbands contain strong contributions 
from highly- i onized iron (the F e XXV/Fe XX complex 
at 0.185 nm; iWhite et a l."2005') as well as free-free con- 
tinuum emission. It is reasonable to assume, therefore, 
that the X-ray emission observed by GOES would em- 
anate from approximately the same flaring plasma for 
which the electron densities were derived. 

An accurate determination of the electron density is 
important for understanding both the heating and cool- 
ing of flare plasmas, and the mechanisms responsible. For 
example, the high densities often inferred from the SXR 
emitting corona during flares are believed to be a result of 
chromospheric evaporation, whereby chromospheric ma- 
terial is ablated up into the corona by a beam of non- 
thermal electrons as inferred from blucshift s of high- 
temperature emission lines /e.g.. [Antonucci fc Dennisl 
IT9831 ICzavkowska et all 119991 . IMilhean fc DennisI 12009^ 
Density values are therefore useful for determining the 
mass rate into (and out of) the overlying loops during 
evaporation (and condensation). Flare cooling is also 
believed to be dominated by thermal conduction around 
the p eak of a flare once th e injected energy is switched off 
(e.g., iRafterv et al. I [20091 ). This transitions to radiative 
cooling at some point during the decay phase. As this 
cooling scales as the square of the electron density inaccu- 
rate assumptions of iVg may lead to large inaccuracies in 
the amount of radiative cooling. Similarly, the amount 
of energy radiated during a flare can be est imated us- 
ing GOES data. Recently, iRvan et all ()2012D measured 
the total radiative losses for ~50,000 flares as a function 
of GOES class. A constant density of 10^° cm^'^ was 
assumed across each event. The results presented here 
show that large flares can have densities 1-2 orders of 
magnitude higher, and that that value can vary during 
individual events. Combining GOES and EVE observa- 
tions can therefore lead to a more accurate determination 
of flare energetics. 

The volumetric fllling factor, /, can also be determined 
by a combination of imaging and spectroscopy once the 
density is known. EM can be derived from continuum 



observations from GOES or the Ram aty High- Energy 
Solar Spectroscopic Imager (RHESSI; |Linet_alJ [l002) , 
while the volume can be estimated from imaging in- 
struments such as RHESSI, GOES/Soft X-ray Imager, 
Hinode/X-ray Telescope or the Atmospheric Imaging As- 
sembly (AIA) also on SDO, although AIA is known to 
saturate during the largest flares. High coronal den- 
sities may also help explain the occurrence of coronal 
HXR emission which is believed to be due to electrons 
accelerated from a coronal reconnection site impinging 
upon the underlying loop. Previous estimates place this 
value as low as 10^-10""^° cn i~^, although these w ere often 
from much weaker events (jKrucker et al.l [20081 ). Precise 
knowledge of the flare loop density can be used to es- 
tablish whether the observed emission is due to thick- or 
thin-target bremsstrahlung. 

Despite being designed to measure changes in the so- 
lar EUV irradiance, several recent studies, in addition to 
the work presented here, have demonstrated how EVE 
observations can be used to determine son ie of the fun- 
damen tal properties of solar flare plasmas: iWoods et al.l 
(j2011h presented evidence f or a flare 'late phase', as well 
as coronal dimming; Hudson et al.l ([2011]) were able to 
derive Doppler velocities despite EV E's modest spec- 
tral resolution; iMilligan et al.l ()2012D showed that it is 
possible to determine the timing and energetics of the 
free-free and fr ee-boun d continuum during large events; 
and rChamberl in et al.l (2012) describes how EVE can 
be used to measure the thermal evolution of the flaring 
plasma. Combining these, and other, analyses of solar 
flares throughout Solar Cycle 24 will give us new insights 
into their global behaviours and properties. 
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